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Abstract. The use of adiabatic passage techniques to mediate particle transport 
through real space, rather than phase space, is becoming an interesting possibility. 
We have investigated the properties of Coherent Tunneling Adiabatic Passage 
(CTAP) with alternating tunneling matrix elements. This coupling scheme, not 
previously considered in the donor in silicon paradigm, provides an interesting 
route to long-range quantum transport. We introduce simplified coupling 
protocols and transient eigenspectra as well as a realistic gate design for this 
transport protocol. Using a pairwise treatment of the tunnel couplings for a 5 
donor device with 30 nm donor spacings, 120nm total chain length, we estimate 
the time scale required for adiabatic operation to be ~ 70 ns, a time well within 
measured electron spin and estimated charge relaxation times for phosphorus 
donors in silicon. 
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1. Introduction 

The potential of silicon based proposals in particular to build upon and integrate 
with existing CMOS fabrication techniques makes them an attractive candidate for 
large scale quantum information processing architectures. Since the original proposal 
by Kane [l] for qubits based on the nuclear spins of '^^P donors in isotopically pure 
silicon ^®Si, other donor based systems have been suggested based on electron charge 
[5] and spin [31 H] [5] degrees of freedom. Scale-up, however, requires more than the 
ability to fabricate many qubits and the associated nanoelectronics. In Ref. fS] a 
scalable donor electron spin qubit architecture was proposed based on a bi-linear 
arrangement of qubits. A key feature of this architecture is a transport mechanism 
using Coherent Tunneling Adiabatic Passage (CTAP) [6] which both allows for qubit 
transport during quantum error correction and serves to reduce the effective gate 
density of the nanoelectronics. CTAP is a spatial analogue of the well known STIRAP 
protocol [7] from quantum optics for coherent transfer of quantum information. The 
advantage of CTAP is that it affords long-range, flexible and robust transfer as is 
required for scale up of donor based architectures. 

In addition to phosphorus in silicon, there are proposals to observe CTAP 
in quantum dots [9l [TOl, superconductors [11], single atoms in optical potentials 
[ladS], and Bose-Einstein Condensates [HIIS- Recently Longhi et al. fTTl IT8] 
demonstrated CTAP of photons in a three-waveguide structure adding significant 
motivation to demonstrations with massive particles. CTAP has also been proposed 
as a means of implementing a quantum analogue to fanout using MRAP (Multi- 
Recipient Adiabatic Passage) [T9l[20j and operator measurements. Other recent works 
have investigated the differences in CTAP in quantum optics and solid state systems 
[211 [22] ■ Atomistic simulations of triple donor systems in silicon have been investigated 
[53] , verifying the existence of an CTAP pathway for a small three donor device and 
the use of ion-implantation for fabrication of CTAP devices have been discussed Ref. 
[24] . Triple dots have recently also been demonstrated in GaAs 2DEG structures 
[25l [26l l27] and carbon nanotubes [28] . 

The transport protocols discussed in [6] [8] were based on the straddling scheme 
introduced by Malinovsky and Tannor [53] where two sites are coupled to the ends 
of a strongly coupled and uncontrolled chain. Providing that the tunneling matrix 
elements between sites on the chain are significantly higher than the gate controlled 
tunneling matrix elements (TMEs) of the end-of-chain sites, then the protocol works 
by transporting the particle between the end-of-chain sites with negligible population 
ever appearing along the intermediate sites in the chain (even transiently). 

In this paper we consider a different adiabatic protocol suitable for long range 
(multi-site) transport based on the alternating scheme. This was proposed in the 
context of quantum dot transport of electrons by Petrosyan and Lambropolous [9] 
and is also known for STIRAP protocols [30] Here we consider this scheme for the 
case of phosphorus ions in silicon. 

Successful coherent transfer of the electron across the chain requires the protocol 
to be completed within the charge relaxation time, and in order to be useful for spin 
based quantum computing architectures must also be within the qubit spin dephasing 
time, T2. Recent measurements of electron spin relaxation for phosphorus donors 
in purified Si^^ have yielded a T2 of approximately 60ms at 7K ^31] and still in the 
millisecond range when the donors are placed near a surface or interface [32j . Charge 
dephasing in a Si:P-P+ system have been investigated in Refs. [33l[34], estimating the 
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Figure 1: (a) Schematic of a five donor, one electron device for CTAP using the 
alternating coupling scheme. Barrier gates (By) control the tunneling of the electron 
between buried phosphorus donor sites, (b) The counter intuitive coupling scheme. 
Tunneling matrix elements in the five site chain are controlled as an A-B chain. The 
counter-intuitive pulse sequence applies all of the fl2 tunneling before fli tunneling to 
effect transport from |1) to |5). 




dephasing time to be of order nanoseconds. 

We consider a five donor structure, illustrated in Figure [U^a) , where four donors 
are ionized and one is neutral. Gates on the surface control the tunnel matrix elements 
and thereby define the adiabatic pathway for electronic transport. The counter- 
intuitive, alternating gate bias sequence, illustrated in figure [TJb), varies the tunnel 
matrix elements resulting in robust transport of the electron from one end of the 
chain to the other. This structure is extendable to an arbitrary (odd) number of sites. 
Based on an effective mass treatment we analyse the system and make estimates of the 
timescale required for the operation of CTAP with the alternating coupling scheme 
for this device. 

2. Analytics 

We begin by writing down the analytical form of the null space for the problem 
with 2n + 1 states for a single particle. We assume that the TME between states is 
\2i — 1) and \2i) is fti, and that all of the even fii and the odd fi^ can be controlled 
separately from to some maximum value. Note that we are here explicitly using 
the modal approximation. Recent calculations have shown this to be a very good 
approximation for CTAP problems [TH [HJ [551 US] although the exact nulling of some 
central state populations that is predicted for CTAP does not precisely occur for 
realistic potentials. Assuming that the energies of the states are all degenerate, we 
write down the Hamiltonian as: 

n 

H = Y^ f^2»-i|2z - l)(2z| + n2^\2^){2i + 1| + h.c. (1) 

2 = 1 

We can immediately write down the null space solution, which we write as |I?o) = 
J2i=i V'iN)- First note that all of the at even sites must vanish, i.e. il^2i — 0. 
Secondly we can write down the values of the ip at odd sites as a series, i.e. 

= - ^-!'~V 2»-1, (2) 
"2i 

and hence 
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We may now write down the normalisation to get the correct form of the null vector, 
|I?o): which is 



1 / " 

n 

+ ■■■+ (-i)"n"2«-ii2n+i 



(-i)^n^2^n^2^-ii2j+i) 

i=j i=l 



(4) 
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N : 



\i=l 



n n 1 

yi=j i=i 



n 



1/2 



(5) 



This form for |Po) immediately shows that this null space has all the desired properties. 
It is unidimensional, when the odd TMEs are all zero and the even TMEs all non- 
zero, the system is in |1), and when the even TMEs are all zero with the odd non-zero, 
the system is in state \2n + 1), with a smooth adiabatic pathway between the states. 
Apart from enforcing the zeros of the tunneling matrix element no control as the 
relative values is required. In fact, this adiabatic pathway is maintained provided that 
only the ends of the alternating pathway can be satisfactorily nulled, ie that fii = 
at t=0 and = at t = tmax, which makes the protocol an attractive prospect for 
implementing as less fine controls are required, decreasing the amount of gates and 
connections required. 

Considering the |I?o) for a five site system (ACTAP5) explicitly we have: 

.„(5), r!2r!4|i)-r!ir!4|3) + r!ir!3|5) 

) = — , =. (6) 

v/(r!2r!4)2 + (f7if74)2 + (r!if}3)2 

We see that when Vli is zero, there is only population in |I?o) Similarly, when Q.^ 

is zero, the null state is |I?o) = |5) does not depend on 0.2- This allows us to consider 
global controls of all the even tunneling matrix elements together and similarly for all 
the odd elements, with only the need to fine tune the ends of the chain, f2i and in 
the ACTAP5 case. Global controls, and specifically A-B chains as we are discussing 
here, have been investigated in the context of Heisenberg chains [35] . 

Another important feature of equation [S] is that it illustrates the robustness of the 
ACTAP protocol to fabrication errors. In particular, one should note the existence 
of the null state irrespective of the values of the TMEs, providing that they can be 
varied between zero and finite values. This is useful when considering the limitations 
of conventional fabrication techniques, e.g. single ion implantation |36j which will have 
unavoidable variability in the final positions of the donors due to straggle. As only 
the ratios of the TMEs enter into the form of the null state, the adiabatic pathway 
persists in the presence of such variations in site to site coupling. It is important to 
note that the adiabaticity of the protocol, and hence the time taken for high-fidelity 
transport, will vary also, and so any given device will need to be calibrated. Further 
discussions on the role of TME variability in related CTAP protocols can be found in 
%2^2Ai. 
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For simplicity let us now consider the case where all the odd tunnel matrix 
elements are the same, and the even ones also, i.e. 

fl2i~i = ill and fl2i = fl2 V«, (7) 

and concentrating on ACTAP5 we have the following result for the null state: 

(5) ^ nl\i)-n,n2\3) + nl\5) 

Furthermore we can also determine all of the eigenstates and eigenvalues for this 
symmetric case. The eigenstates are 

(5) ^ »i|l)-»il]2|3)+r!f|5) 

\vf) = T Jni-nin2 + ni{\2) - 14)) 



+ (f^i-02)|3)+r!2|5)][2.^r!2-r!i02 + r!i]-\ (10) 
\vfl) = [f^ill) ± ^^1 + ^1^2 + ^1^ + |4)) 



+ (r!i + r!2)|3)+r!2|5)][2^r!2 + r!i02 + r!2]-\ (11) 

with eigenenergies 

=0, (12) 



E± ^ ± J VLl - VLin2 + VLl, (13) 



E2± = ± ^^1 + ^1^2 + 9.1 (14) 

To explore the evolution of the system, for analytical simplicity we chose square 
sinusoidal pulse variation of fii and $72, however in keeping with most adiabatic 
protocols, the exact form of the variation in tunneling matrix elements is not essential. 
For pulses from < = to i = <„iax we set: 



= f^max sin^ ( 7^7— ) ) (15) 

^^2 = r!„iaxC0s2 (-^) , (16) 

and these pulses are illustrated in Figure [^^a). 

The null state is the eigenstate which has zero energy throughout the protocol. 
This is the state used for the CTAP protocol, as its evolution under the pulse sequence 
gives rise to a smooth change in population from |1) at t = to |5) at i = tmax- We can 
then find the eigenvalues for this Hamiltonian, shown plotted in figure [2l^b) . Solving 
the master equation for this Hamiltonian, we obtain the occupancies of each of the 
donor sites throughout the protocol. Analogously with other CTAP protocols, we 
treat the electron as initially occupying the first site, and transferred entirely to the 
end-of-chain site at the end of the protocol. In the alternating scheme there is transient 
occupation of the site in the middle of the chain, as shown in Figure Wis) but other 
intermediate sites remain unoccupied. This differs from the straddling scheme [6l [8] 
where occupation of all intermediate sites is strongly suppressed. 
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Figure 2: (a) The top trace shows the TMEs fli and O2 as a function of time. Below, 
the TMEs between donors (postions marked by dashed hues) are illustrated as a 
density plot as a function of time. TMEs controlled by gates 2 and 4 are high initially 
(^2), whilst the TMEs controlled by gates 1 and 3 are low initially (r^i). Varying 
the biases on the gates effects the counter-intuitive pulse sequence described, (b) 
Eigenvalue spectrum of the ACTAP5 Hamiltonian . (c) Populations of the position 
eigenstates throughout the protocol. Note the complete transfer from |1) to |5) 
but with transient population in |3), and no population in cither |2) or |4). (d) 
Adiabaticity as <max increases. As imax is increased the adiabaticity parameter is 
decreased indicating better fidelity tranfer across the chain. A shorter total time for 
the protocol can to transitions out of the jPo) state and the transfer of the electron 
from |1) to |5) is no longer complete. 



The Adiabaticity parameter can be fairly easily calculated between the null state 
and either of the nearest neighbours, so in particular, choosing the closest positive 
state, we write the Adiabaticity parameter as 

^_\<p^mMi (17) 

^ r2i»i-»2»f + »1^2(»l-»2) , . 

For adiabatic evolution we require that ^ ^ 1 with the maximal adiabaticty 
determining the time allowable for a realistic experiment. If we note that the 
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adiabaticity parameter will be largest at the crossing point of the two TMEs, and 
if we apply pulses that are symmetric, then we can make a considerable simplification 
to the above. So at the point where A is greatest, we have fi^ = = r2max/2, 
and rii = D,2 = 7''f^max/(2tmax) and we find using the form of the TMEs as given in 
equation I15[ the above equation gives rise to the simple form 

A^-^ . (19) 

V Oi imax^max 

This dependence of the adiabaticity parameter on the the total protocol time ^maxis 
plotted in Figure [D^d). With longer tmax, and hence lower A, the transported electron 
is more likely to remain in the desired \'Dq) state resulting in better fidelity transfer. 

A realistic device may not have enough fine control over the TMEs to be able to 
completely turn them on or off. In this case we have non-zero minimum values and 
assuming a form for the TMEs of: 

f^i (t) = f7i „iin + {^1 max ^ ^1 min) siu^ ( j , (20) 

^2{t) = f72min + (f^2max " f^2min)c0S^ ( oT" ) ' (^^) 

\ ^^max / 

In the case where complete suppression of the tunneling is not possible, a measure 
of our ability to turn off tunneling between donors is the contrast ratio between the 
overlap of the initial and final states with |I?o)- The contrast ratio defines an error 
rate for the protocol when we do not have complete suppression of the tunneling. 
For a device to successfully perform the protocol we require the contrast ratio to be 
high, reaching a maximum of 1 when complete suppression is possible. This can be 
determined by calculating the product of the overlap (II'Dq) at i = and (I'olS) at 

t — ^max- 

(1|-no(/ = n)) = _ ^ "2max .32) 

V 1 mill ' 2 max ' 1 mm 2 max 

and assuming r^imin ^ ^2 max, to first order: 

(l|I?o(i = 0)) = (23) 

h max 

Similarly, at the end of the protocol, when t — imax the null state is: 

(VM = <,nax)|5) = Imax .34) 

V 2 mm ' Imax ' 2 mm 1 max 

Applying the same approximations as before gives the result: 

(Po(i = W)|5) = l-;^s^, (25) 

^"l max 

SO a first approximation to the fidelity of population transfer in this case of imperfect 
contrast ratio is 



|2 _ 1 ^lmin^2 

002 o2 

max^'2max 



|(5|I?o(t = W))(I?o(t = 0)|l)r - 1 - ZTW"""^^- (26) 



To illustrate with an example, to achieve a fidelity of 99.9% and assuming that 
f2i = ri2, we require a ratio of rimin/^^max = 0.3, hence we can see that a high 
fidelity transfer may still be possible, even with imperfect controls. 
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Figure 3: a) An interdigitated gate design for A-CTAP. Donors are placed in between 
gate "fingers" in the central region of the device. All odd gates are connected together 
and similarly for all the even gates, b) A slice taken from a TCAD simulation of 
potentials within the device, slicing through the centre of all the donors which are 
located in a line parallel to the position axis. 



We now turn our attention to calculating the gate modulated tunnel matrix 
elements for a semi-realistic 5 donor device which we will use with the above 
expressions for adiabaticity and fidelity to estimate the timescale of the protocol over 
such a device. 

3. Numerical Modeling 

The increased simplicity of an alternating coupling scheme over other schemes makes 
it an attractive candidate for possible implementation for demonstrating CTAP in a 
Si:P setting. Figure [3] illustrates a possible gate design to implement ACTAP5. This 
design features global control of the even and odd tunneling matrix elements with 
even and odd barrier gates connected to be pulsed together. These global controls 
are desirable due to the simpler fabrication and reduced control circuitry required. 
Fine tuning of the end of chain sites may be implemented using further gates but they 
are not included here. Modeling of the potentials within the device with applied gate 
biases was performed using a commercial TCAD finite element modeling package. A 
simulation of the correct gate voltage sequence required to enact the CTAP protocol 
is a difficult control problem with cross-talk effects needing to be taken into account. 
Compensation schemes for gated donor schemes have been investigated in Ref. [37] 
and similar measures may need to be implemented here to ensure the correct tunnel 
matrix elements are obtained. We apply here a simpler approach of looking at pairwise 
tunneling rates with a simplified structure consisting of two donors and a barrier gate 
to modulate tunneling rates. 

The phosphorus in silicon scheme we are considering consists of ionised ^^P 
buried in a ^^Si substrate. Each of the ionised phosphorus atoms forms bonds with 
neighbouring silicon atoms with their remaining 4 valence electrons. One un-ionised 
donor remains in the chain with a remaining valence electron loosely bound to the 
donor. This remaining electron's wavefunction may be manipulated using surface 
gates to effect the transport across the chain. 
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Figure 4: (a) Ungated As^s as a function of donor separation in [100]. As donors 
are moved further apart the tunnehng rate between them is decreased. The large 
fluctations at large donor separations are due to fiuctations in the Monte Carlo routine, 
(b) A slice taken from a TCAD simulation of the simplified, one gate, two donor device. 
The slice in the vertical plane of the two donors. The 10 nm wide gate is centred about 
250 nm directly between the positions of the two donors placed 30 nm apart, as marked 
by the crosses in the picture. (Inset) Relationship between fitted charge density and 
TCAD obtained fields. Note the linear fit. Errors in the value of A obtained were 
±0.01%. 



The electronic states of shallow donors in silicon have been well studied. To 
calculate the energy states of the remaining valence electron we use the effective 
mass approach developed by Kohn and Luttinger [38] in which the donor electron is 
described using hydrogenic envelope functions. This approach is extended to include 
a second ionised donor such that the electron is described by symmetric and anti- 
symmetric superpositions of the ground state for each donor and an external electric 
field. The effective Hamiltonian for the two donor, SiiP-P"*" system is given by: 

n = nsi + Vd{r - + Vd{r - R2) + Ve (27) 

where iJg; is the Hamiltonian of an electron in the pure silicon lattice, which includes 
both a kinetic term and the effective potential due to the silicon lattice, Vd the 
coulombic potentials of each donor and Ve an external applied field [39l [40l [41] . The 
valence electron can be described by symmetric and anti-symmetric superpositions of 
the ground state of each for each separate donor: 

^-i (r) = M[4^{v - R„) ± i^(r - R^)] (28) 

where A/" is a normalisation factor and the the ground state for a single phosphorus 
donor centred about Ra is given by: 

6 

^(r) = F^{r - i?)e*'=''-(''--«°). (29) 

M = l 

The F^{r — R) are effective mass, hydrogenic envelope functions about the six 
degenerate band minima. 

TMEs may then be calculated by examining the minimum energy splitting 
between the lowest energy states of a donor pair, the symmetric-antisymmetric gap 
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^SAS ■ TMEs may then be varied from fimax = ^SAS to using time varying gate bias 
pulses such as those described in equations [13 and [111 Examining these gate driven 
coupUngs gives an indication of the timescale of the tunnehng rates, and hence the 
timescale of the protocol. 

Figure [4] shows the calculated Asas as a function of donor spacing for a P-P+ 
system with no external field is applied. At 30nm Asas is calculated as O.lfieV, 
corresponding to a tunneling matrix element of ~ 2 GHz. These calculations are 
consistent with the hydrogenic results of Openov [42] To include the effect of the 
external field on the tunneling rate we have utilised an analytic potential fitted to 
the potential landscapes obtained from several TCAD simulations at different gate 
biases. This approach allows us to use a simple form for the external field which 
varies smoothly over space which can be varied with any possible gate potential while 
maintaining a close connection to modeling of semi-realistic structures. In the y-z 
plane we approximate the surface gate with an infinite line of charge the same distance 
above the ground plane with the origin at the line of charge. Using the method of 
images we obtain the form of the potential: 

T/^-^in ^;+;-+j;; m 

27reoer + (z - dy 

where A is the charge density of the line of charge obtained using the simulated TCAD 
potentials, shown in figure [4l[b)Inset and d is the distance of the line of charge above 
the ground plane. Since the dielectric constant of the Si02 is much less than that of 
the bulk silicon, the discontinuity in the resulting field is accounted for by extending 
the 5nm oxide layer to an equivalent thickness with a bulk silicon dielectric constant to 
ensure the correct field in the bulk silicon is obtained. This relationship between A and 
barrier gate voltage has a constant offset due to the contact potential of around 0.5 V 
that arises from the aluminium gates placed on the silicon dioxide layer. Since we 
are only interested in the effect net effect of the gate we have ommitted this constant 
offset in the calculation of Asas so zero volts corresponds to the field-free (ie purely 
hydrogenic) case. In practice this would be achieved by the gate offsetting the surface 
charge. 

Figure [5| shows the calculated variation of Asas with barrier gate voltage. We 
see a variation in Asas from ~ 0.5 to 10 GHz with a change of approximately 250 mV 
applied voltage on the barrier gate. This is consistent, within uncertainty, with 
previous modelling of the energy gap for P-P+ charge qubit with donor spacing of 
30nm, showing that tunneling may be varied from around to 10 GHz with a change 
in barrier gate bias of 300 mV . The errors in Asas dominated by fluctuations 
in Vz due to sampling errors in the Monte Carlo routine. Using equation [^Sj with 
^min — 0.5GHz and fimax = lOGHz and assuming f^i = ^2 we obtain an error rate of 
10^*^. This result demonstrates the remarkable robustness of the protocol to imperfect 
control of the tunneling matrix elements. Ideally Vz is when the field is completely 
symmetric with respect to the two donors and consists of the diagonal elements of the 
external field component of the Hamiltonian cancelling with each other. If we neglect 
Vz in calculating Asas we obtain a minimum tunneling rate of around 0.1 GHz and an 
error rate of 10~^. 

Assuming perfect contrast, for a target adiabaticity of 0.01 we may now calculate 
the required time for the protocol. Using equation I19i with r2inax=10GHz we obtain 
imax ~ 70ns as the required time for adiabatic transfer of an electron across the chain, 
in this case 120nm long. 
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Figure 5: (a)Variation of Agas witli applied barrier gate voltage. We see a variation in 
Asas from 0.5 to 10 GHz in a change of approximately 250 mV applied voltage on the 
barrier gate. The larger errors at larger spacings are due to the Monte Carlo routine, 
(b) Vx and Vz coefficients as a function of barrier gate voltage. The fluctations in Vz 
are due to sampling errors in the Monte Carlo routine. 

Allowing for a large range variablity in the TMEs while still being able to obtain 
high fidelity transfer of the electron means that the protocol is also robust to donor 
placement. Variations to donor placement will result in variations in the ungated 
TMEs as shown in Figure ID However, as shown in equation [261 adiabatic behaviour 
will still occur over a wide range of maximum and minimum TMEs. As a result 
the protocol can be insensitive to variations in donor placement, which is important 
since devices such as this proposed are generally fabrication using ion implantation 
techniques where there is variability in the exact placement of the donor. 

4. Conclusion 

We have shown that for a semi-realistic device consisting of buried phosphorus donors 
in silicon with surface gates, an adiabatic pathway may be found allowing transport of 
quantum information along a donor chain using the alternating coupling scheme. The 
fidelity of the transport protocol is resilient to variations in the TMEs along the chain 
so long as the end-of-chain sites can be controlled. We calculate that, for a device with 
a relatively large donor spacing of 30nm, with a maximum pairwise tunneling rate of 
10 GHz which can be suppressed to 0, we can estimate a time of ~ 70 ns required for 
adiabatic operation. For closer positioned donors we expect an adiabatic operation 
time down to the few nanosecond regime or faster. This time is well within measured 
electron spin relaxation times and dephasing for silicon and within the expected charge 
relaxation times for Si:P devices. 
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